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ABSTRACT 

The observed strong dark-to-luminous matter coupling is described by a bias relation 
between visible and dark matter sources. We discuss the bias which emerges in the 
case where the topological structure of the Universe at very large distances does not 
match properly that of the Friedman space. With the use of such "topological" bias, 
we construct the Universal Rotation Curve (URC) for spirals which occurs to be in a 
striking agreement with the empirically known URC. We also show that the topological 
bias explains the origin of the Tully-Fisher relation (L ~ V 4 ) and predicts peculiar 
oscillations in the URC with a characteristic length ~ \JL. 

Key words: galaxies: kinematics and dynamics - galaxies: spiral - dark matter. 



1 INTRODUCTION 

It has been long known (Persic, Salucci & Stel 1996, PSS 
hereafter) that the shape of rotation curves of spirals is 
rigidly determined by a single global parameter, e.g., lu- 
minosity or the number of baryons in a galaxy. This feature 
was stressed in PSS by an empirical construction of a uni- 
versal rotation curve (URC) which describes quite well the 
rotation velocity at any radius and for any galaxy as a func- 
tion of, say, the galaxy luminosity only. It follows that the 
distribution of the Dark Matter (DM) in galaxies carries a 
universal character as well, and is a function of the luminous 
mass. Note that the standard Cold DM (CDM) models fail 
to explain this strong dark-to-luminous matter coupling, by 
an obvious reason: in any model where DM is built from 
hypothetical non-baryonic particles (e.g., CDM, worm DM, 
or self- interacting DM) the number of the DM particles in 
the halo of a galaxy is, essentially, a free parameter and re- 
lating it to the number of baryons in the galaxy requires 
some very strong nonlinearity. Moreover, it is well estab- 
lished that the DM density in galaxies shows an inner core, 
i.e. a central constant density region (e.g., see Gentile et al., 
2004 and references therein; Weldrake, de Blok & Walter 
2003; de Blok & Bosma 2002; for spirals and Gerhard et al., 
2001; Borriello, Salucci & Danese 2003 for ellipticals), which 
is in a clear conflict with the predictions of A -CDM models 
yielding NFW-type profiles with a cusp (e.g. pr>M ~ 1/r) 
in the central region of a galaxy (Navarro, Frenk & White 
1996). 

The strong coupling between DM halos and baryons 
(see also Donato, Gentile & Salucci 2004) definitely requires 
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some new physics. The coupling can be described by a rigid 
relation between the sources of dark, pdm, and visible, pl, 
matter, the so-called bias relation (Kirillov 2006). In the 
linear case the most general form of the bias relation is 



Pdm (x) = b (x,x) pl (x) dV'. 



(1) 



The homogeneity assumption b (x, x') — b (x — x') allows 
one to fix empirically the bias operator b em p- Indeed, in this 
case the Fourier transform of the bias relation Q gives 



Pdm (i, k) — b (i, k) p L (i, k) 



(2) 



where we added a dependence on time to account for 
the cosmic evolution. The empirical bias function b emp = 
Pdm (t,k) /pl {t,k), in virtue merely of its definition, will 
perfectly describe DM effects at very large scales (i.e. in the 
region of linear perturbations) . 

The present Universe is not quite homogeneous though, 
e.g. it is not uniform at galaxy scales. Still we would expect 
relation J3J to hold in the geometrical optics limit (i.e. for 
rather short wavelengths as compared to the Hubble scale, 
or to a cluster scale when a single galaxy is considered). 
Parameters of the bias function may then vary for differ- 
ent spatial regions, i.e. b emp may include an additional slow 



dependence on the location in space: b : 



, (t, k,x). In or- 



der to fit observations, any theoretical source of DM should 
reproduce properties of the bias function b emp in detail. 

In the linear gravity, the bias relation @ can be inter- 
preted as a modification of the Newton law: 



! + &(*)) 



sin (kr) — krcos(kr) 
kr 2 



dk. 



(3) 



The asymptotically flat rotation curves in galaxies require 
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that the correction to the Newton's potential should be log- 
arithmic, i.e. the gravitational acceleration should switch 
from r~ to r _1 . This, according to 10, implies b(k) ~ fc _1 , 
or 



b(x — x ) ~ \x — x |" 



(4) 



at galaxy scales. In fact, observations suggest the same be- 
havior of b(x — x') for much larger scales (Kirillov 2006). 
Indeed, the distribution of the luminous mass shows char- 
acteristically fractal behavior: the mass Ml(t) within the 
ball of radius r grows, essentially, as ~ r D with D ~ 2 on 
distances up to at least 200Mpc (Pietronero 1987; Ruffini, 
Song & Taraglio 1988, Labini, Montuori & Pietronero 1998) . 
The bias function ||1J leads then to Mdm{t) ~ r 3 , i.e. to a 
homogeneous distribution of the DM and, hence, of the to- 
tal mass. Thus, bias @ is just one that reconciles the two 
seemingly contradictory observational facts: the fractal dis- 
tribution of baryons with the dimension D ~ 2 and the 
large-scale homogeneity of the metric. 

A theoretical scheme capable of explaining the origin of 
such bias was proposed in (Kirillov 1999, Kirillov & Turaev 
2002). It was shown there that processes involving topology 
changes during the quantum stage of the evolution of the 
Universe unavoidably (and, in fact, model-independently) 
lead to a scale-dependent renormalization of the constant of 
gravity G, and this effect can be imitated by the emergence 
of DM, whose distribution is linearly related to the distri- 
bution of actual matter. Importantly, assuming the thermal 
equilibrium during the quantum stage predicts in almost 
unique way a very specific form of the bias function (Kir- 
illov & Turaev 2002; Kirillov 2003) 



b(k) = 



/' 



Vn 2 + k 2 



for k < (j,. 



(5) 



where fi ~ Tpia (tpi) /a (t) has the meaning of the primor- 
dial temperature at which the topology has been tempered 
and k ~ ma (tpi) /a (t) is the mass of primordial scalar par- 
ticles if they do exist. This means that at scales k < k the 
bias becomes the constant b (k) = jx/k and the standard 
Newton's law restores. However at galaxy scales it should 
be k ~ jj, 2> K and in what follows we take k = in O 
1 . Thus in the coordinate representation the bias takes the 
form 



Hr,t) = — 7 / {b(k)k 



2tt 2 



, 3n sin (kr) dk 



kr k 



(6) 



/' 



27r 2 r- 



(1 — cos (fir)) . 



Bias J5J is of form 0, so it predicts the logarithmic cor- 
rection to the Newton's potential for a point source: 8(f) ~ 
j£- lnr at r > Ro where Ro — ir/ (2/x) (see for details Kir- 
illov & Turaev 2002; Kirillov 2006). Thus, the parameter R 
plays the role of the scale at which DM starts to show up, 
so in galaxies it has to be estimated as a few Kpc. As bias 
10 has a thermodynamical origin, there have to be certain 



1 E.g., inflationary scenarios suggest m ~ 10 — 5 mp; and, there- 
fore, fi/n ~ 10 5 . However so far scalar particles have not been 
observed which makes us think that scalar fields are not more 
but phenomenological objects. 



fluctuations in the value of Ro (this effect is analyzed in the 
next Section). 

In the present Letter we demonstrate that bias 0, JSJ 
gives a very good agreement with the empirical URC con- 
structed in PSS and, therefore, it seems to be reasonable 
to believe that the nontrivial topological structure of the 
Universe gives indeed a correct explanation to the DM phe- 
nomenon. 



2 THE TULLY-FISHER RELATION 

The parameter /i in JHJl has the meaning of the tempera- 
ture at which the topology has been tempered (Kirillov & 
Turaev 2002; Kirillov, 2003). We expect that in the very 
early Universe it had the order of the Planck temperature 
Tpi (recall that at T > Tpi quantum gravity effects are 
thought to dominate). When, on the further stages of the 
Universe evolution, topology changes were suppressed, this 
parameter evolved as fi ~ Tpia(tpi) /a(t), where a is the 
scale factor. Note, however, that the present value of /j, has 
the sense of the primordial temperature. As it is extremely 
small: fi ~ 10 _23 T 7 where T 7 is the temperature of CMB 
radiation (T 7 ~ 2.7 K), we have to admit the existence of 
a specific phase in the past when the nontrivial topological 
structure might decay (Kirillov 2002, 2003), causing a cer- 
tain re-heating of matter. During the decay phase, \ia was a 
decreasing function of time, i.e. the scale Ro = n/ (2/i), that 
corresponds to the cross-over from the standard Newton's 
law to the logarithmic behavior of the potential of a point 
mass, grew faster than the scale factor a(t). 

Note that the homogeneity of the Universe requires the 
total mass distribution (luminous plus dark components) 
to have a constant density in space. With the bias of form 
0, @, this corresponds to a fractal distribution of baryons 
(Kirillov 2002, 2003), i.e. the number of baryons within the 
sphere of a radius R > Ro behaves as 



N b (R) ~ vR L 



(7) 



with D ~ 2, while for R < Ro the fractal distribution is 
unstable (for the Newton's law restores and baryons dom- 
inate over the DM). The increase of Ro/a allows one to 
assume that in the very early Universe there was a moment 
i« when Nb (Ro (£*)) < 1, i.e. baryons had the fractal dis- 
tribution ||7J at all scales. After the topology decay phase, 
as the scale Ro (t) "jumps" towards a new, higher value, the 
fractal distribution is preserved at scales larger than Ro, but 
it becomes unstable on smaller scales. The instability devel- 
ops and baryons under a certain scale of order Ro start to 
redistribute, governed by Newtonian dynamics. This means 
that we can relate Ro to the scale of galaxy formation. Then, 
according to Q, we should expect the number of baryons 
in a galaxy to be 



N b ~ u'BS 



(8) 



with the values of Ro and v' corresponding to the moment 
when a galaxy started to form. Note, however, that during 
the formation of a galaxy the value of Ro switches off from 
the Hubble expansion (Kirillov 2006), i.e. law JHJ remains 
valid for the present-time values of Nb and Ro. 

It is easy to see that this relation leads directly to the 
Tully-Fisher law L ~ V 4 (Tully & Fisher 1977), where 
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L is the luminosity and V is the rotation velocity of a 
galaxy. Indeed, recall that _Ro fluctuates in space: small 
spatial fluctuations of the primordial temperature p (A/U//U 
~ AT 7 /T 7 ) represent seeds for the present-time scatter in 
the local value of fj, — 7r/ (2Ro) in different galaxies. Accord- 
ingly, the masses of galaxies M ~ rribNi, (where rut is the 
baryon mass) fluctuate as 



M ~ mtf Ro ■ 

This fixes the choice of 



2 \ m^' 



-l/D 



(9) 



(10) 



in the bias JSJ , © for any given galaxy. Sufficiently far from 
the center, the galaxy can be considered as a point-like ob- 
ject, so ©,© yield the following law for the gravitational 
acceleration at a sufficient distance from the edge of the op- 
tical disc: 



(11) 



GM 9 (l , 2 / • I \\\ 

9 = —pr- [ 1 + n ^ r ~ sin ^o r >> ) ■ 

As V^a/r = g{r), the Tully- Fisher relation for the asymp- 
totic rotation velocity Vac follows now from l|ll|l . ill Oil : 



Vl ~ -GM a u 



L a 



M a 



V^ 



(12) 



with (3 = J2_ and a 2 = §G (m b u') 1/D . Thus, in our inter- 
pretation of the DM, the Tully-Fisher law reduces to rela- 
tion (1 1 ( )l) which, in turn, can be read as an indication of the 
fractality of the primordial distribution of baryons with the 
dimension D ~ 2 (sec Q). 



3 ROTATION CURVE OF SPIRALS 

Let us now compute the rotation curve (RC) of a galaxy 
modeled by an infinitely thin disk with surface mass density 
distribution pl — ae^ r ' R °5{z). From Q, ||S|), we find for 
the DM halo density (we use the notations x — r/Rn and 
A = /iRd) 



Ph (x) = 



Act 



(1 — cos (A \x — y\)) d y (13) 



2tv 2 R d J \x-y\ 2 

where y lies on the plane z = 0, while x is the 3-dimensional 
vector. For the sake of convenience we present the Fourier 
transform 

M L 



If we neglect the oscillating term and compare this with the 

2 

pseudo-isothermal halo p = po -j", i we find for the core 
radius 



2 _ Re 



R 2 D In (1 + A 2 ) ■ 



(17) 



According to PSS, the core radius can be estimated as 

a = 4.8(L/L») 1/5 (18) 

with log L» = 10.4, which makes A a certain function of the 
luminosity. 

Consider now circular velocities predicted by the above 
mass distributions. For the disk contribution to the equilib- 
rium circular velocity, we get PSS 



V 

— ¥ = f D {x, A) = 

J2 



(19) 



- T- ("(f) *(!)- Ml) *(?)) 



and for the Dark Halo contribution we find from I14H the 
expression 



= f H (x,X) = - 



o (yiPTiy) 



j Ji (kx) dk, 



(20) 



where J„, I n , K n are the Bessel and the modified Bessel 
functions and V 2 = D a -A. Thus for the rotation curve 

00 Rjj n 

we find the expression 

V 2 (x, A) = Vl (f D (x, A) + f H (x, A)) . (21) 

As we see, the shape of the rotation curve depends indeed 
on one parameter A. Via relation J5J, or equivalently 112L A 
is expressed as a function of the total number of baryons 
Nb in the galaxy; there is, however, an uncertainty in A 
due to the variation of the ratio M/L for different galaxies. 
At the moment of the galaxy formation Rd ~ -Ro which 
corresponds to the same initial value A ~ 1 in all galaxies. 
On subsequent stages of the evolution A becomes different in 
different objects. Indeed in smaller galaxies supernovae are 
more efficient in removing the gas from the central (stellar 
forming) region of a galaxy than in bigger galaxies (e.g., 
see Shankar et al., 2006 and references therein) and this 
creates the fact that in smaller objects the disc has a smaller 
baryonic density (a lower surface brightness) and the ratio 
A ~ R D /R > 1. 

To compare expression 1211 with that from PSS we 
rewrite it as 



Ph (k,k z ) = 



/' 



Vk 2 + ki (( fcflx ,)3 + i)i 



{jx - x/PTfef) (14) V 2 (x,\) _ f D (x,\) + f H (x,\) 



* opt 



f D (3.2, A) + f B (3.2, A) 



(22) 



where 6 is the step function: 6 (u) — for u < 0, and 9 (u) = 
1 for u > 0, and Ml = 2-k<jR d is the (non-dark) mass of 
the galaxy. 

First of all we note that this distribution is quite con- 
sistent with the observed cored distribution (Gentile et al., 
2004). Indeed, in the central region of the galaxy 

Ml 



Ph(0) 



(2tt) 2 R D 



Aln 1 + A" 



while for i>lwe find 

2p H {0) (1- cos (Ax)) 

pH(a) " ln(l + A 2 ) x 2 ' 
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(15) 



(16) 



where x = 3.2 corresponds to the optical radius of a galaxy. 
The plot of this curve for different values of A is presented 
in Fig.l. 

While the similarity of our RC 12211 with the empiri- 
cal URC of PSS is quite good, we note that the topological 
bias JSJ predicts a new feature in Rotation Curves - spe- 
cific oscillations in the DM density with the characteristic 
wavelength t = 2-7r/A (or in dimensional units t ~ M g ). 
Indeed at a sufficient distance from the edge of the optical 
disc (i.e., as x S> 3.2) a galaxy can be considered as a point- 
like object. Then from lilt for the rotation velocity we find 
the expression 
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Figure 1. The rotation curves V 2 (x) /V 2 t vs x = r/R]j for different values of A . The green and red dashed lines give the visible and 
DM contributions respectively, while the black line gives the sum. We see that with the decrease of the luminosity (increase of A) DM 
fraction increases in agreement with PSS. 






K 

2Xx 



sin (Ax) 
A.r 



(23) 



which shows the presence of a specific oscillations with the 
decaying amplitude l/(Xx). Such oscillations are, in turn, 
rather difficult (though possible) to extract from observa- 
tions. Indeed in the case of HSB (high surface brightness) 
galaxies when A < 1 (i.e., for rather long periods £ > 2n) 
the expression fl'-j.'ili gives a very good quantitative approx- 
imation to the exact formula 12211 starting already from 
x — x opt — 3.2. However the reliable RC data available 
extend usually not more than to x — (2 — 3)x op t- In this 
range oscillations are not established yet and the beginning 
of oscillations is seen (e.g., see Fig.l.) as not flat RC slopes. 
The slopes observed are known to take values between 0.2 
and —0.2 (e.g., see PSS). In the case of LSB galaxies A ^> 1 
(£ <C 27r), the amplitude of oscillations is somewhat sup- 
pressed ~ 1/A and the small amount of the stellar mass 
in the range x op t < x < 3x opt considerably smooths such 
oscillations which results in a some deviation of (123ft from 
the exact expression 12211 . Moreover, in deriving 1221 we do 
not take into account the presence of gas which due to su- 
pernovae does not trace the brightness, i.e., it deviates the 



exponential profile. Essentially this is true for LSB galaxies. 
Thus, to observe such oscillations we have either to measure 
velocities for sufficiently large distances ~ lOxopt (e.g., for 
HSB), or to improve the accuracy of available observational 
data in LSB. 



4 DISCUSSION AND CONCLUSIONS 

As it can be seen from Fig.l, the topological bias JSJ pre- 
dicted in (Kirillov & Turaev 2002; Kirillov 2003) shows quite 
a good agreement with observations. Indeed, it repeats all 
features of the empirical URC of PSS: the amount of DM 
progressively increases with the decrease of the luminos- 
ity (cf. PSS), DM shows the cored distribution (cf. Gentile 
et al., 2004) with the strong correlation 117H between the 
core radius and the disk size (cf. Donato, et al., 2004), and 
the Tully-Fisher relation (Tully & Fisher 1977) is explicitly 
present (see 1121 1. There is no doubt that the tuning of a 
single free parameter M/L allows to fit any RC. At least, it 
is claimed in (e.g., see Milgrom & Sanders 2005 and refer- 
ences therein) for the RCs obtained via Milgrom algorithm 
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MOND (Milgrom 1983), and our RCs are phenomenologi- 
cally quite close to those, although the physics in our ap- 
proach is completely different. In this respect we can claim 
that the topological bias gives a rigorous basis for applying 
the MOND-type algorithm in galaxies (which however al- 
lows for the Tully-Fisher relation to have j3 ^ 4). Therefore, 
there is enough evidence that bias Q , (JSJ gives an adequate 
description of galaxies. 

We repeat that our approach produces as good fit to 
the observed RCs as the Milgrom algorithm (known to be 
quite successful empirically (Milgrom & Sanders 2005, see 
however Gentile et al. 2004, Donato et al. 2004)) can do. 
However, contrary to MOND, our theory remains linear in 
weak fields, and the superposition of forces holds. In fact, 
our approach does not presume any modification of the the- 
ory and basic equations: there is actually no modification 
of gravity, while the bias appears merely as a result of a 
disagreement between the actual topology of the physical 
space and that of the flat space (e.g., see Sec. 2 in Kirillov 
2006). Thus, there are all reasons to believe that the DM 
phenomenon has indeed the topological origin. 

Once we accept the topological origin of the bias Q, 
(|SJ, the Tully-Fisher relation 



V r 



8 



(24) 

with /3 = -§^i — 4 serves as a strong indication of the frac- 
tal behavior in the primordial distribution of baryons with 
the dimension D ~ 2. Such fractal distribution changes es- 
sentially the estimate for the baryon number density in the 
Universe (e.g., see Kirillov 2006). The currently accepted 
post-WMAP cosmology has (roughly): Qtotai — 1, fiA ~ 0.7, 
fioM ~ 0.25, and flj ~ 0.05, which implies QoAi/^b ~ 5. 
We stress that such estimates are model dependent, for they 
are strongly based on the standard model (e.g., the con- 
tent, evolution, the homogeneity of the baryon distribution, 
the power law of initial spectrum of perturbations, etc.) 2 . 
Moreover, the direct count of the number of baryons gives 
fi;, ~ 0.003 for the whole nearby Universe out to the ra- 
dius ~ 300/i^g 1 Mpc (e.g., see Persic & Salucci, 1992), which 
means that in the standard cosmological models most of 
baryons are somewhere hidden. 

When the topological bias is accepted such estimates 
require essential revision. Indeed, according to © the topo- 
logical bias modifies the Newton's law in the range of scales 
p > k > k where the equilibrium distribution of baryons has 
the fractal behavior. On scales k < K the standard Newton's 
law restores (and baryons cross over to the homogeneity) but 
dynamically every particle becomes heavier in 1 + p/n times, 
which gives for the effective DM fraction fioM /fit ~ 1+p/n. 
The scale Ro = ir/ (2/x) is directly measured in galaxies by 
RCs and is estimated as a few Kpc. However the mean value 
< Ro > for the homogeneous Universe should be 10 2 times 
bigger (Kirillov 2006). The maximal scale l/« is the scale 
where the primordial fractal distribution of baryons crosses 
over to the homogeneity. This scale is not so easy to mea- 
sure without a detailed investigation. Indeed, the large-scale 



2 For linear perturbations there are no doubts that if we take the 
observed spectrum of AT/T and fix an arbitrary model a(t), we 
find in a unique way the initial primordial spectrum of Ap/p. To 
fix the model we have to know both the measured AT/T and the 
initial values of Ap/p which a priori are unknown. 



structure, e.g., the existence of huge (~ 100 — 200Mpc) voids 
with no galaxies inside and thin (~ 1 — 5Mpc) walls filled 
with galaxies, fits quite well into the fractal picture and sug- 
gests only the lower boundary l/« > 100 — 200A4pc. This 
gives for the DM fraction fioM /fib > 10 2 — 10 3 which is 
consistent with the observed value fij, ~ 0.003. However the 
maximal possible value 1/ft ~ Rh (Rh is the Hubble ra- 
dius) which gives fiuivi/fit, ~ 10 cannot be excluded. To 
avoid misunderstanding we stress that the topological na- 
ture of the bias makes the fractal distribution to be equi- 
librium and consistent with the homogeneity of the metric 
and the observed CMB fluctuations AT/T (e.g., see Kir- 
illov 2006). The topological nature means that the same 
bias appears in all interactions. If the bias would not mod- 
ify the electromagnetic field, then the fractal distribution 
of baryons would be in a severe conflict with observations 
and surely had to be rejected as it does take place in the 
standard models (e.g., the fractal distribution produces too 
strong fluctuations AT/T ~ Apb/pb ~ /x/ft ). The topo- 
logical nature of the bias however creates the fact that the 
Coulomb force and all Green functions are also modified at 
galaxy scales (Kirillov & Turaev 2002, Kirillov 2006) which 
reduces AT/T to the observed value Ap to tai / ptotai ~ k/ ' p 
~ 10 -5 (e.g., for sufficiently remote objects < p > r ^> 1 
the apparent luminosity has to behave as I ~ L/r 15-1 in- 
stead of 1/r 2 which gives for the number of objects brighter 
than £ somewhat higher (with respect to D — 3) value 



N (1) w vr D (I) « 1/e 



D/(D-1) 



). Thus the topological bias 



and the observational definition of 1/k requires the careful 
and thorough revision of the standard model and all basic 
formulas. 

In conclusion, we point out that bias (|SJ predicts the 
existence of specific oscillations in the distribution of DM 
with the characteristic wavelength ~ M g ~ y/TJ. When 
the observational data allow, this can be used to verify the 
theory and, thus, to make a more definite conclusion on the 
nature of DM. 
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